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Q> ■ As a study of peculiar velocities of nonlinear structure, we analyze the model of 

' a relativistic thin-shell void in the expanding universe. (1) Adopting McVittie (MV) 

0^ ■ spacetime as a background universe, we investigate the dynamics of an uncompensated 

void with negative MV mass. Although the motion itself is quite different from that of a 
^p^i compensated void, as shown by Haines & Harris (1993), the present peculiar velocities 

O I are not affected by MV mass. (2) We discuss how precisely the formula in the linear 

I perturbation theory applies to nonlinear relativistic voids, using the results in (1) as well 

' as the previous results for the homogeneous background (Sakai, Maeda, & Sato 1993). 

(3) We re-examine the effect of the cosmic microwave background radiation. Contrary 
^ to the results of Pim & Lake (1986, 1988), we find that the effect is negligible. We 

show that their results are due to inappropriate initial conditions. Our results (l)-(3) 
suggest that the formula in the linear perturbation theory is approximately valid even 
for nonlinear voids. 
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1. Introduction 

Measurements of large-scale peculiar velocities can provide a constraint on the universe model 
(see, e.g. Zehavi & Dekel (1999)). In contrast to other cosmological tests, they give a constraint 
on the density parameter alone, almost independent of the cosmological constant Aq (Carroll, 
Press, & Turner 1992). 

Although the relation between the peculiar velocity and the density parameter J7o is usually 
given in the linear perturbation theory (LPT) (Peebles 1976), the observed universe has nonlinear 
density profiles. In fact, a network of nonlinear voids filling the entire universe has been suggested 
by redshift surveys such as the CfA2 (Geller & Huchra 1989) and the SSRS2 (da Cost et al. 1994). 
Moreover, using such redshift surveys, the description of a void-filling universe was confirmed 
quantitatively (El- Ad & Piran 1997; El- Ad, Piran, & da Costa 1996, 1997). The relation between 
and peculiar velocities inside underdense regions suggests Qo < 0.3 can be ruled out at the 
2.4-sigma level (Dekel & Rees 1994). 

It is therefore important to investigate peculiar velocities of nonlinear void structure. Here 
we consider the model of a relativistic thin-shell void. The expansion law of relativistic voids 
was investigated originally by Maeda & Sato (1983a,b), developing the metric junction method 
proposed by Israel (1966). They found analytically that, in the flat universe, the shell radius R 
expands asymptotically as i? oc t(i5+Vi7)/24 ^0.797 (Maeda & Sato 1983a). For the self-similar 
void model, on the other hand, Bertschinger (1985) obtained the solution with R oc The 
difference of the two results is so small that we do not have to care which model is better, which 
is physically determined by the radiative process in the shell. For other universe models, the 
motion of the shell was calculated numerically (Maeda &: Sato 1983b): in the open universe, the 
shell expansion is eventually frozen to the background expansion; on the other hand, in the closed 
universe, the shell expands much faster and its velocity finally approaches the speed of light. The 
relation between the peculiar velocity of the shell and the universe model was later investigated 
systematically (Sakai, Maeda, & Sato 1993). 

Lake & Pim extended the work of Maeda & Sato so as to include a mass inside a void (Lake & 
Pim 1985) and the cosmic microwave background (CMB) radiation (Pim &z Lake 1986, 1988). In 
particular, they claimed that the inclusion of radiation has significant quantitative and qualitative 
effects on the evolution of the void. It was shown, for instance, that the asymptotic behavior of 
the shell is i? oc t in the flat universe if the CMB radiation is included (Pim & Lake 1986). This 
result is in contrast to that for the vacuum void (R oc t^'^), and hence quite surprising. 

Haines & Harris (1993), on the other hand, included a mass outside a void by employing 
McVittic (1966, hereafter MV) metric instead of Friedmann-Robertson- Walker (FRW) metric. MV 
metric approximates a spherical mass embedded in an asymptotically FRW spacetime. "MV mass" 
represents the degree to which the void is not compensated by the mass of the shell. They demon- 
strated the history of the shell in the flat MV spacetime, showing that the negative MV mass acts 
to accelerate the shell expansion. 
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In this paper we extend the previous work to clarify the following points. 

(1) Peculiar velocities of uncompensated void, which is characterized by negative MV mass. 
Although Haines & Harris (1993) discussed the dynamics of voids with non-zero MV mass, the 
effect of MV mass on the peculiar velocity is not clear. It is important to see it because there 
is no evidence that actual voids have the shells with compensated mass. For example, if voids 
originate from primordial bubbles that are nucleated in a phase transition during inflation (see, e.g. 
Amendola et al. (1996)), it is unlikely that voids have compensated shells. 

(2) The relation between and peculiar velocities was derived in LPT (Peebles 1976). It is 
important to see how precisely the formula applies to nonlinear relativistic voids. We address this 
question, using the results obtained in (1) as well as the previous results (Sakai, Maeda, & Sato 
1993). 

(3) As we mentioned above, Pim & Lake (1986, 1988) arrived at the surprising conclusion 
that the effect of CMB radiation is significant. If it is true, we should take it into account seriously 
when we constrain the cosmological parameters from the observation of bulk flow. Therefore, their 
result deserves closer examination. 

This paper is organized as follows. In section 2, we present the relativistic equations of motion 
for a thin shell in MV spacetime, which will be solved later numerically. In section 3, we investigate 
how the peculiar velocity changes due to MV mass for uncompensated voids. In section 4, we 
compare the results for the present model and those in the LPT. In section 5, we examine the effect 
of the CMB radiation on the void evolution. These results are summarized in section 6. 

2. Basic Equations 
2.1. McVittie Spacetime 

The spacetime described by MV metric has several useful properties: 

• The near-field limit is Schwarzschild, in isotropic coordinates; 

• The far-field is a FRW spacetime; 

• The energy-momentum tensor has a perfect-fluid form. 

It is an exact embedding of the Schwarzschild metric into the FRW metric. 
The line element is 

ds' = - (\^) ' dt^ + a\t){l + h)\dx' + f{x){d0' + sin^ Od^^)}, (2-1) 
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where p and p are the energy density and the pressure, respectively, which are inhomogeneous in 
general. H = {da/dt)/a is the Hubble parameter at % — oo- m is a constant and called MV mass. 
The h ^ limit is clearly FRW, while Schwarzschild solution is recovered by a — ^ 1. If A; = or 
— 1, the scale factor a{t) always takes the Priedmann solution. 

To see the meaning of MV mass, let us calculate the local gravitational mass defined by Misner 
k Sharp (1964): 



M=^{l-gi'''d^Rd,R) with R = 



For each k we find 



47r Q 

M{R) = —pR^ + my{x) with y(x) = 



cos5(x/2) {k = +1) 
1 (A; = 0) 

cosh5(x/2) {k = -l) 



(2-5) 



(2-6) 



As long as the void's size is much smaller than the horizon scale, x ^ 1 ^^id y(x) ~ 1 for any 
background model. In this limit, we may therefore interpret the MV mass as approximately the 
Misner-Sharp mass minus the background mass (47r/3)pi?^. 

For a thorough discussion of the McVittie metric, its history, and its place among inhomoge- 
neous models, see Krasinski (1997). 



2.2. Junction Conditions 

Let us derive the equations of motion for a spherical shell around a void, by developing the 
thin-shell formalism of Israel (1966). The basic equations for the shell in the flat MV spacetime 
were given by Haines & Harris (1993). Here we rewrite the equations as in a simpler form, which 
also describe the shell in a closed or open background, by extending the equations of Sakai, Maeda, 
& Sato (1993). Because we are interested only in the effect of the outer MV mass, we assume the 
inside region to be homogeneous throughout the paper. 

Let a time-like hypersurface S, which denotes the world-hypersurface of a spherical shell, divide 
a spacetime into two regions, (outside) and V" (inside). We define a unit space- like vector AT^, 
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which is orthogonal to S and pointing from V~ to V~^. It is convenient to introduce a Gaussian 
normal coordinate system (n,x') ^ in such a way that the hypersurface of n = corresponds to 
S. From the assumption that a shell is infinitely thin, the surface energy- momentum tensor of the 
shell is defined as 



5^1, = lim 



a [ T^^dn. (2-7) 



Using the extrinsic curvature tensor of the hypersurface of the shell, Kij = N^j, and the 
Einstein equations, we can write down the jump conditions on the shell as (Berezin, Kuzmin, & 
Tkachev 1987) 

[Kij\^ = -SttG (^Sij - ^hijTrS) , (2-8) 
= Kf, (2-9) 

' ^ ' Si = (2-10) 

where hij denotes the three metric of S, and | denotes the covariant derivative with respect to hij. 
We have denoted the value of any field variable on S on the side of by and defined a 
bracket as [*]^ = — ^~ . Eliminating K-~ from equations (2-8) and (2-10), we can derive the 
equation (Berezin, Kuzmin, & Tkachev 1987): 

Kj^Si + ^G{s^Si - liTvSf} = [T«]±. (2-11) 

If the outer region is homogeneous, equation (2-11) leads to a simple expression of the basic 
equation, because it does not contain the metric in V~ Sakai & Maeda (1993a,b); Sakai, Maeda, 
&: Sato (1993). On the other hand, if we assume the inside region to be homogeneous, which is 
the case here, it is easier to solve the equation obtained by eliminating K-j~: 

KfSf - 47:G{s^Si - \{TTSf} = [T„"]±. (2-12) 
The line elements in y+ and in V~ are described by 

/ -, _ T \ 2 

= - [iTh) '^*+ ^ + h)Hdxl + flix+)id0^ + sin^ ^V)}, (2-13) 

ds'^ = -dtl + a'i{t^){dx'i + fl{x-){de^ + sin2 ^V), }, (2-14) 
The direct calculation of yields (Sakai &; Maeda 1993a) 

„_ 7_(/' +v-H_R) ^ ndv- , , 



^In this paper, Greek letters run from to 3, while Latin letters run in 0, 2, 3. 
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where the circumference radius of the shell R, the peculiar velocity of the shell v- , and its Lorentz 
factor 7_ are defined as 

dx- 1 
R = = a-f-, V- = a--r—, and 7_ = , (2-16) 

Similarly, u+ and 7+ in the MV spacetime are defined as 

_ a+{l + hf dx+ , 1 
^+ = — ] 7 — -TT- and 7+ = (2-17) 

As energy-momentum tensors, we consider perfect fluid on S and in V^, i.e., 

Sf,u = (cr + ct)i;^?;;;^ + zuh^i,, (2-18) 

= iP^ + P^Hu^u^ + (2-19) 

where a, vj, v^, and arc the surface density, the surface pressure, the four velocity of the shell, 
and the four velocity of the background fluid, respectively. In the Gaussian normal coordinate 
system, we have 

r-±=72(t;2p + p)|±, T^^=j^v{p + p)\^. (2-20) 

Now, with the help of equations (2-15), (2-18), and (2-20), we can write down equations (2-9) 
and (2-12) explicitly as 

da dR a + vo , , , , 



= -7-|(l-2-J.-iJ--^-|-2.0(. + 4^)- -' - ; . (2-22) 
The relation between dR/dt- and v- is given by 

^ = f'_v_ + H_R. (2-23) 
Further, the conditions of the continuity of the metric, 

dT^ = (YTjTfc) *5- - »+('+)(! + = <it- - «-(*-)<JX-, (2-24) 

f = A{(i + ,)V/,} = ;|:(a_/_). (2-25) 

reduce to 

^ = ^, (2-26) 

dt- l-h-i- ^ ' 
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The equations of motion for the shell are determined by equations (2-21), (2-22), and (2-23). We 
use the above supplementary equations (2-26) and (2-27) to give t+ and f+, respectively; they are 
used at the initial time as well as at each step of time evolution. 

The angular component of the jump condition (2-8), 

7+ (/; + + v+H+R^ - 7- ifL + v-H^R) = -AirGaR, (2-28) 

gives a constraint for the relation between the surface density a and MV mass m. We use it for 
giving initial data as well as for checking numerical errors of integration. For integration we adopt 
the 4th order Runge-Kutta method. Throughout the analysis wc did not encounter any numerical 
problem: the relative errors of equation (2-28) were always less than 10^"*^'^. 

The equations of motion presented here and by Sakai, Maeda, & Sato (1993) have several 
advantages compared with the equations derived in other papers. First, the expression is much 
simpler. Secondly, there is no sign ambiguity in the relation between t+ and t-, (2-26), contrary 
to the comment by Pirn & Lake (1986). Thirdly, our expression for the extrinsic curvature Kg, 
(2-15), can take both positive and negative values without ambiguity. Although our equations make 
numerical integration easier, they may not be so convenient for analytic arguments. 



3. Peculiar Velocities of Uncompensated Voids 

Here we consider only dust as matter fluid: 

= and w = 0, (3-1) 

where the subscript oo denotes quantities at x+ ~^ oo. A compensated void simply means m = 0, 
i.e., the background is described by the FRW metric. On the other hand, an uncompensated void 
is characterized by negative MV mass. Here we fix the value of m by supposing no shell (cTj = 0) 
at the initial time tj. The initial time ti is a free parameter, which is determined by the structure 
formation model; in the following we set ti as the decoupling time , i.e., Zi = 1000. The remaining 
initial parameters are fixed as follows: 

v+ = 0, R,H+ = 0.1, H+ = H-, pr=0, 

= = 1 or 0.98. (3-2) 



Figure 1(a) shows the motion of the shell in terms of the comoving coordinate x- As shown by 
Haines & Harris (1993), negative MV mass pushes the shell faster. Although Figure 1(a) indicates 
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that the effect of MV mass looks quite large, the behavior of % is not observable. What we can 
observe is the radius and velocity of the shell at the present epoch. We thus plot the peculiar 
velocity normalized by the Hubble expansion: 



in Figure 1(b). The asymptotic behavior is determined by Qq, independent of MV mass. Figure 
2 reports the relation between vq and Qq, which confirms that vq does not depend on whether the 
void is compensated or not. 



4. Comparison with Linear Perturbation Theory 

In this section we discuss the relation between Qq and vq, comparing the results in the rela- 
tivistic void model and those in LPT. 



The peculiar velocity v for general density fluctuations in LPT is (Peebles 1976) 

with FRiQ'^-^, (4-1) 



3HQ 

where g is the peculiar gravitational acceleration. For a spherically symmetric system, the gravi- 
tational acceleration is given by 

9{R) = (4-2) 

where 6M{R) is the difference between the mass within a sphere and the unperturbed mass within 
the sphere with the same radius R. 

For the void model, 5M{R) depends on whether we measure it just inside the shell {R = R-) 
or just outside it {R = R+): 

Air Air 
SM{R^) = -^{p+ - p-)R\ SM{R+) = -^{p+ - p-)R^ + AiraR^ (4-3) 

It is therefore reasonable to define the mass difference as the average: 



On the other hand, one of the junction conditions (2-28) can be rewritten as 



(4-4) 



= -A^GcjR, (4-5) 

where = signi^^^. In the Newtonian approximation, {dR/dr)^ <C 1, y{x) ~ 1) and = +1, 
equation (4-5) reduces to mass conservation: 

m + — (/9+ - p-)R^ = AnaR^. (4-6) 
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Prom equations. (3-3), (4-1), (4-2), (4-4), and (4-6), we obtain 



6 V P+ 47rp+i?V3y 
For compensated voids (m = 0), equation (4-7) reduces to a simple expression: 



1-^1. (4-8) 



6 V P 

For uncompensated voids (cTj = 0), on the other hand, equations (4-7), (4-6) and p'^a\ =const. 
read 




(4-9) 

Figure 3 shows plots of vq v.s. J7o for compensated voids, where the subscript denotes 
quantities at the present. (The details of the analysis for the homogeneous background were given 
by Sakai, Maeda, & Sato (1993). In the linear case (a), our numerical result is in good accordance 
with the result in LPT. Even in the nonlinear case (b), the difference between the two results is 
relatively small (up to 10%). 

Let us turn to the case of uncompensated voids. Obviously the term [rf /r^)^ in equation 
(4-9) represents the effect of MV mass: as the comoving radius r increases, the effect of MV mass 
decreases. This argument explains the result that the eventual behavior of v does not depend on 
MV mass, as shown in Figure 1(b). 



5. Effect of CMB radiation within a void 

As we mentioned in the introduction, Pim &; Lake (1986,1988) showed that, if we include 
CMB radiation, the shell expands much faster than that in the absence of radiation, and that its 
asymptotic behavior is i? oc t even in the flat universe. Here we re-examine the effect of radiation 
in the flat FRW background: m = and = 0. 

First, let us reproduce the results of Pim & Lake (1986). As background matter, a mixture of 
dust (p^) and blackbody radiation (p+) is considered: 

p^=p^+pt=Pcr^^^, p-' = 'f (5-1) 

Setting the present temperature Hubble parameter as Tq = 2.7K and Hq 
using the relation, 

+ _ Stt^ k%T^ 
~ 15 h^c^ ' 



= 100/kms/Mpc and the 



(5-2) 
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the background model is completely fixed. The interior is assumed to be the flat FRW spacetime 
with radiation only, of which abundance {p~) is characterized by a parameter, 



For matter fluid on the shell, they assume that the equation of state has a form, w = ea. In one 
of their calculations, the initial parameters arc flxcd as follows: 



Integrating the equations of motion in section 2, we obtain the result in Figure 4(a), which is the 
reproduction of Figure 5 of Pim & Lake (1986). 

This flgure tells us that, no matter how little radiation exists, it affects the shell's motion 
significantly. Because this result was surprising to us, we examine their analysis. As a result, we 
find that their assumption of k~ = was inappropriate for the following reason. If fc^ = and 
> , the Friedmann equation reads > H~ . As Sato (1982) and Sato & Maeda (1983) 
argued, however, a thin shell is formed by compression of matter like a snow-plow mechanism 
when the inner expansion is faster than the background expansion, i.e., H~ < . Therefore, the 
assumption of A;^ = is inconsistent with the thin shell description. If we still used the thin-shell 
equations for the case where the shell expands faster than the interior matter fiuid, a part of the 
shell mass would be forced to "evaporate" so as to keep the inside homogeneous, i.e., the shell 
would emit mass and accelerate, which seems unphysical. This explains the odd behavior in Figure 



Here we re-analyze voids with CMB radiation. Although the exact value of /H^ cannot 
be determined without knowing the formation process, the consistency with thin shell requires 
H~ > . Because we still assume the background universe to be flat, the inside region should be 
open. Adopting — instead of = and leaving the other conditions unchanged, we solve 
the equation of motion. The result is reported in Figure 4(b), showing that the effect of radiation 
is much smaller. 

We should note, however, that it is not so fruitful to investigate further details of the motion of 
voids including radiation in this approach. Because we do not know the physical process of radiation 
around the shell, the equation of state (e) is not determined; furthermore, even the validity of the 
thin-shell approximation is not clear. What we can conclude is that, under the condition that the 
thin-shell approximation is valid, the effect of CMB radiation on void expansion is negligible. 




(5-3) 



Zi = 1000, vl = 0.1, R,H+ = 0.1, e = 0, k 
a = 10-\ 10"^ 10-^ 10"^, or 10"^ 



= 0, 



(5-4) 



4(a). 



6. Summary 



As a model of nonlinear structure, we have considered a relativistic void in the expanding 
universe, and discussed peculiar velocities. 
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(1) In order to investigate the dynamics of shells with uncompensated mass, we have adopted 

McVittie spacetime as a background universe. Although the motion itself is quite different from 
that of a compensated void, as shown by Haines & Harris (1993), the present peculiar velocities 
are unaffected by MV mass. 

(2) We discuss the relation between Qq and peculiar velocities, comparing the results in the 
present model with those in the linear perturbation theory. For nonlinear voids, the quantitative 
difference between these two results is up to 10%, which is relatively small. 

(3) Because Pim & Lake (1986,1988) arrived at the surprising conclusion that the effect of a 
small amount of CMB radiation is significant, we have re-examined it. We have shown that their 
results are due to inappropriate initial conditions. With modified initial conditions, the effect of 
radiation turns out to be negligible. 

Although we have investigated only specific models of nonlinear structure, our results (l)-(3), 
as a whole, indicate that the formula for peculiar velocities in the linear perturbation theory can 
apply approximately to nonlinear voids. 

Numerical Computation of this work was carried out at the Yukawa Institute Computer Fa- 
cility. N.S. was supported by JSPS Research Fellowships for Young Scientists, No.9702603. 
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FIG. 2. Relation between vq and for compensated (m = 0) and uncompensated (m < 0) 
voids. 
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FIG. 4. Effect of CMB radiation, (a) Reproduction of Fig. 5 of Pirn & Lake (1986). Both 
the outside and inside are assumed to be flat FRW. Results for various a are presented, (b) The 
outside is flat, but the inside is open: = , a = 0.1. 
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